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ABSTRACT 



We present HerschelfPA.CS, MMT/Hectospec and XMM-Newton observations of Abell 1835, one of the brightest X-ray clusters on 
the sky, and the host of a strong cool core. Even though Abell 1835 has a prototypically "relaxed" X-ray morphology and no signs 
of ongoing merger activity in strong- and weak-lensing mass maps, it has a complex velocity distribution, suggesting that it is still 
accreting significant amounts of mass in the form of smaller satellite systems. Indeed, we find strong dynamical segregation of star- 
forming dusty galaxies from the optically selected cluster population. Most Herschel sources are found close to the virial radius of the 
cluster, and almost a third appear to be embedded within a filament feeding the cluster from the SW. We find that the most luminous 
infrared galaxies are likely involved in galaxy-galaxy interactions that may have triggered the current phase of star formation. 
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1. Introduction 

Galaxy clusters have long been thought of as particularly ac- 
tive sites of galaxy evolution (iDressleJ 11980 ). Early optical 
studies suggested that clusters affect stellar populations in in- 
falling substructure mainly by turning off star formation as 
galax ies encounter the dense cluster environment (Lewis et id] 
' l2002h . However, when these same objects were observed in 

■ the infrared by IRAS, ISO and Spitzer, a significant fraction 
of cluster galaxies was reveale d to contain h e avil y obscured 
sites of star formation (see e.g. iMetcalfe et al.l (1200 5) and ref- 
erences therein). When this obscured activity is accounted for, 

. clusters emerge as environments where star formation is both 

■ quenched and induced, perhaps even at different stages of the 
same physical process. Some studies have claimed to observe 

■ trigg ered star form ation in strongly disrupted, merging clusters 
' (e.g. iGeac h etaP 12006) but the intrinsic scatter in this correla- 
tion is large, with some seemingly relaxed clusters presenting 
much higher star formatio n activity than their actively merging 
counterparts (Hai nes et al.l f2009). 

Many studies in recent years have attempted to separate clus- 
ter samples into two distinct populations: merging clusters vs 
clusters that are dynamically relaxed. The global properties of 
the latter are expected to be tightly coiTelated with mass, and 
therefore directly relatable to statistical predictions from cosmo- 
logical models. 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



While there is evidence that some X-ray cluster proper- 
ties may in fact be bimoda l, e.g. with respect to the presence 
or absence of a cool core dSanderson et al.l l2009h . cosmolog- 
ical simulations have shown that even the most relaxed cool 
core clusters can s till be accreting significant amounts of mass 
(iPoole et al. 2008), at a mean rate which is app roximately lin- 
early proportional to the mass of the system dMcBride et al.l 
l2009h . Generally, only the most massive of these accretion 
events (mass ratios < 10 : 1) are expected to disturb the dense 
cluster core. Massive cool core clusters can thus quietly accrete a 
significant fraction of their total mass without large disturbances 
to their central ICM, which makes them arguably "cleaner" tar- 
gets for studies of cluster galaxy evolution. 

The Local Cluster Substructure Survey (LoCuSSQ) Herschel 
key programme (.Smith.2010.) was designed to probe dusty star- 
forming galaxies in the infall regions of a sample of 30 galaxy 
clusters (z ~ 0.2) spanning a wide range of mass and merging 
histories. 

Abell 1835 is the most X-ray luminous cluster in th is sample 
(Lx.boi = 5.32 + 0.15 X 10'*^ ergs/s (Zhang et alj|2008|)) and the 
host of a strong cool core (Peterson et al. 2001). Its strong- and 
weak-lensing derived mass distributions (Richard et al. 201 
lOkabe et allBOlO.) . and X-ray morphology (S mith et aL. 20051 
are typical of undisturbed (i.e. non-merging) cool core clus- 
ters. For example, the low substructure mass fraction in its core 
(/sub = 0.13 + 0.01) sugg ests that it has grown i n mass by < 10% 
in the previous few Gyr dSmith & Tavloiil2008l) . 



http://www.srbham.ac.uk/locuss 
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In this Letter, we focus on a sample of 100/zm-selected 
spectroscopically confirmed cluster galaxies in Abell 1835. We 
use the Photode tector Array Camera and Spectrometer (PACS) 
(|Poghtschl2010l) on the Herschel satellite (Pilbratt 2010) to mea- 
sure far infrared fluxes and from these we estimate total infrared 
luminosities for these objects. We then use these data to examine 
the spatial and velocity distribution of obscured star formation 
sites in Abell 1835. We assume //(> - 70km s"', Q^ - 0.3, D.\ - 
0.7. All virial masses and radii are stat ed with respect to the v irial 
overdensity, AviriaiCz = 0.25) = 121 dBrvan & Normanlll998h . 

2. Observational data 

PACS data were obtained in scan map mode on December 24th 
2009, at both 100 and 160yL(m. The details of the observing mode, 
data reduction p rocedure and source extraction and photometry 
can be found in (ISmithl2010t) . The 90% completeness flux limits 
are 21mJy at 100/im (86 sources) and 28 mJys at 160;um (91 
sources). 

We conducted a wide-field spectroscopic survey of this field 
with Hectospec at MMT in February 2009, confirming clus- 
ter membership of 439 gal a xies ( Hardegree-Ullman et al., in 
prep.). In addition, ICzoskd (|2004|) obtained redshifts for 630 
cluster members with VIMOS/VLT, bringing the total number 
of unique redshifts for this cluster to 851. This extensive spec- 
troscopic coverage means that, out of 91 Herschel detections 
with fidO/jm > 28mJy, 64 have redshifts, of which 33 lie within 
6v < 6000 km/s of the cluster redshift. The total number of spec- 
troscopically confirmed cluster galaxies with > 3cr detections in 
both lOOjjm and 160/im bands is 45, and we plot their location 
in Fig. 1. 

X-ray maps can be used both to separate the AGNs from the 
starburst galaxies and to identify galaxy groups around clusters, 
and we have analysed the XMM observation ID 0147330201 
to characterize the X-ray appearance of A 1835. We decom- 
pose the 0.5-2 keV combined XMM image into unresolved 
and extended emissio n, following the wavelet technique of 
iFino^uenov et al.l (2009), before extracting point source fluxes 
at the locations of Herschel sources. Additionally, we make 
use of previously published wide-field data from Spitzer fMlPS 
24fim and UKIRTAVFCAM J,K bandsl'H aines et al.ll2009i) . and 
Subaru/Suprime-CAM (.Qkabe et al...20101 

3. Results 

3.1. Infrared luminosities, temperatures and SFRs 

At the redshift of the cluster, z - 0.25, the PACS 100 and 160yum 
bands are expected to straddle the peak of the spectral energy 
distribution for sources with dust temperatures of T ~ 30K. The 
160;um/100;um color should thus be sensitive to dust temperature 
differences, with redder objects having cooler dust components. 
We plot observed 160;um/100;um vs 100//m/24yL(m colors for the 
sources in our cluster in Fig. 2, and overplot the syn thetic colors 



from commo nly used templates f rom the li terature: iRieke et af 
120091 m9), 'Dale & Helou"2002 (DH02), iLagache et al.l 12003 



(L03),|Chary & Elbaz 2001 (CEOl). All the SED templates are 
poor representations of these galaxies, which appear too red in 
their 100yL(m/24//m colors. The CEOl templates appear to fit the 
colors of ~ 1/2 of the sample, but with templates corresponding 
to total infrared luminosities that are one to two orders of mag- 
nitude fain ter than the ac tual measured luminosities - see Smith 
( 1201 Ol) and lRawlg ( 1201 Ol) for a more extended discussion of this 
lOOyum excess. 
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Fig. 1. Location of Herschel detected, spectroscopically confirmed, 
cluster members on tlie sky, overlayed on a K-band luminosity den- 
sity image of the cluster produced from all spectroscopically confirmed 
members and corrected for the spatial incompleteness of our spectro- 
scopic survey. The black solid contours correspond to 4cr wavelet de- 
tections of X-ray extended sources in the XMM data. The dashed circle 
marks Tyir = 2.93 Mpc, the dotted square the limits of the PACS field. 
The sizes of the points scale linearly with Ltir, and their color scales 
with the temperature of the best fit modified blackbody SED, ranging 
from 20K for the coldest galaxies to ~ 60K for the hottest. Galaxies 
identified as AGN through emission lines or X-ray emission are marked 
with a green square. 



We estimate total infrared luminosities, Ltir, for these galax- 
ies in two ways: first, we fit the CEOl templates (redshifted 
to the cluster rest- frame) to the 24, 100 and 160 yum emission 
and integrate the best-fit spectrum, normalized to the observed 
fluxes, from rest-frame 3-1100/im. In order to obtain estimates 
for the dust temperatures in these systems, we also fit a mod- 
ified blackbody spectrum to the lOOyum and 160//m emission: 
/v oc [1 - exp(-Tv)]Bv, where By is the Planck function, and 
Ty is the frequency-dependent optical depth, v/v(f. We assume 
yS = 1.5 as the emissivity power law index, and vq = 3 THz 
as the frequency at which the emission becomes optically thick, 
but note that the integrated luminosity depends very weakly on 
the values assumed for /3 e [0.7 - 2] and vq, while the best fit 
temperature varies at most by 20%. The Ltir derived from the 
blackbody fits are systematically lower than the ones obtained 
from the best-fit templates by a factor of ~ 0.6 - 0.8, since they 
underestimate the mid-IR contribution to the total IR flux. The 
cluster members range in Ltir from 4 - 79 X lO'^L©, and all but 
one have dust temperatures in the range 20 - 45K. 

The BCG is much brighter than the rest of the galaxies 
in the cluster, with fluxes of 255.6 + 3.9 mJy (160yL(m) and 
297.8+3.1 mJy (lOOyum) and an estimated Ltir of 7.9 x lO", 
consiste nt with previous Spitzer/MIPS observations of Abell 
1835 (E gami et al.! '2006). We derive star formation rates from 
the combined Hq, (from our spectroscopic me asurements) and 
Ltir emission using the following relation from jKennicutt et al] 
(1^9,) : SFR[Kroupa](M0yr-i) = 5.5 x lO-'^^[L(Ha)obs + 
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Fig. 3. Radial profile of K-band (black) and total infrared (red) lumi- 
nosity surface density relative to the average surface density in concen- 
tric annuli within r <3 Mpc (excluding the BCG). The dashed vertical 
line marks rvi,. The K-band luminosity density is an approximate tracer 
of stellar mass, and declines continuously with radius, whereas the in- 
frared luminosity peaks close to the virial radius of the cluster. 



Fig. 2. 160yum/100/jm vs 100/jm/24yum colors, compared with synthetic 
colors derived from SED templates in the literature. The lOO/jm/24/jm 
colors of the cluster galaxies appear significantly reddened compared to 
the templates. The BCG (yellow triangle) has the lowest 100//m/24//m 
color of the sample, and is well fit by a Dale & Helou starburst SED with 
ax 1 .65. Error bars are based on statistical noise, and do not include the 
overall PACS calibration error. The stated PACS calibration uncertainty 
of < 10% at 100/im is nevertheless not large enough to reconcile the 
data with the template colors. 



0.002 Ltir] (ergss"'). This yields SFRs in the range: 2-20 
Moyr'^ for the general cluster population. The BCG's high lu- 
minosity and special location at the center of a massive cool core 
separate it from the rest of the cluster population and we there- 
fore exclude it from the subsequent analysis. 

We find that, of all galaxies with far-infrared emission de- 
tected by Herschel, only 2 have X-ray counterparts down to 
Lx ~ 4 X 10"^' ergs s"', which is below the classical 10"^^ ergs 
s ' threshold for AGN definition. Only four galaxies (including 
the couple of X-ray point sources) in the sample are spectro- 
scopically classified AGN according to their position in the clas- 
sic BPT diagram (Baldwin et al. 1981), defined by the emission 
line flux ratios of [Nil] 6587/H„ versus [OIII] 5007/H^. We thus 
conclude that the bulk of the infrared luminous population is be- 
ing heated by star formation and not an active nucleus, and we 
exclude the 4 likely AGN from our analysis. 



3.2. Probing the environments of obscured star formation 

It is apparent from Fig. 1 that, apart from the BCG, there are few 
luminous infrared galaxies in the central regions of the cluster. In 
fact, the radial surface density profiles in Fig. 3 show that most 
of the infrared luminosity is concentrated at large (r ~ 2.5Mpc) 
radial distances, which is not the case for the K-band luminosity 
density which, as an approximate tracer of stellar mass, declines 
smoothly with radius. 

Upon closer inspection, however, it is clear that this over- 
density at the virial radius does not occur isotropic ally around 
the cluster, but rather peaks along specific directions. The veloc- 
ity distribution of K-band selected cluster galaxies (black line 




Fig. 4. Redshift probability density function for galaxy populations ob- 
served at different wavelengths. The density function is weighted by 
flux/luminosity, and corrected for the incompleteness in our spectro- 
scopic survey. K-band: black, MIPS 24/jm:blue, Ltir: red 



in Fig. 4) is distinctly non-Gaussian, and shows signs of sig- 
nificant substructure. The most luminous infrared sources, lo- 
cated toward the south of the cluster, are dynamically segregated 
from the rest of the cluster population, clustering around an elon- 
gated structure in the foreground (6v ~ 2000 km/s) of the cluster, 
which is also traced out in the K-band luminosity density maps 
(see Fig. 1). A gradient is observed in the velocities projected 
along the axis of elongation, where the relative velocity with re- 
spect to the cluster center progressively decreases as the inner 
region of the cluster is approached. 

Figure 5 shows Subaru /'-band cutouts of the most luminous 
infrared galaxies, in decreasing order of Ltir- Many galaxies 
appear to have undergone recent interactions — half of them 
have close companions, and half present distorted morpholo- 
gies indicative of recent merging. We can attempt to quantify 
the importance of mergers in our small sample by comparing 
the fraction, /^ , of Herschel sources with close (Ar < 50 kpc 
and Av < 300 km/s) companions to that of a population of non- 
IR luminous cluster galaxies that has been matched, both in K- 
band luminosity and radial distance from the cluster center, to 
the IR luminous sample (excluding in both samples all identi- 
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The Ltir density profile for Abell 1835 is not spherically 
symmetric: the brightest galaxies are clustered together in an 
elongated structure towards the south of the cluster that is dy- 
namically segregated from the rest of the cluster The structure's 
K-band morphology and velocity gradient toward the cluster 
core suggests that these galaxies may be embedded within a fila- 
ment that is still currently feeding the cluster with new, actively 
star forming galaxies. 

XMM imaging of this field reveals the presence of two 
groups in the outskirts of Abell 1835 (see contours in Fig. 1), 
one of which appears to be embedded in the filament that is also 
feeding the bright LIRGs into the cluster. The group is domi- 
nated by a lai-ge (Lk » 2 x 10'^Lo,z = 0.245) elliptical galaxy 
with no detectable Herschel emission. 

The two groups have X-ray luminosities of 4.6 + 0.6 and 
7.5 + 0.8 xlO"*^ ergs s"' (0.1-2.4 keV), from which we estimate 
m asses, using the weak lensing calibration of the Lx-M relation 
of iLeauthaud et al.1 (l20io.) . of M200 = 4.3 and 5.8 x IO'^Mq re- 
spectively. Together these amount to > 1/1 0th of the mass of the 
central cluster, reinforcing the view that Abell 1835 continues to 
grow by feeding on smaller systems from the surrounding large 
scale structure. 
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Fig. 5. Subaru i' band ~ 60 x 60 kpc image cutouts centered on the 
IR brightest cluster galaxies, in decreasing order of logjg ^tir (upper 
right comer). A significant fraction has close companions (c), and half 
show low surface brightness tidal features, or asymmetric morphologies 
indicative of a recent interaction (D) 



fied AGNs). We find that, while the total sample of IR lumi- 
nous galaxies is indistinguishable from the matched sample of 
IR faint galaxies (/,. ~ 12+^% vs 10+2%, respectively), the IR 
bright (logjoLxiR > 10.8) galaxies in Abell 1835 are marginally 
more likely to have a close neighbour (/c ~ 30 
the respective matched sample). 
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4. Summary and discussion 

Star formation can be triggered by several diff'erent physical pro- 
cesses prevalent in clusters and their infall regions; disentangling 
these eff'ects is one of the main aims of the LoCuSS Herschel key 
programme. The most promising candidates, ICM ram pressure, 
galaxy mergers, and galaxy harassment can all perturb the ISM 
to produce a new generation of stars, but the relative importance 
of these proc esses i s a stro ng function of clustercentric radius 
(c.f. Fig. 1 in lSmithI (l2010h ). The location of the peak of IR lu- 
minosity density just inside the virial radius of Abell 1835 sug- 
gests that merging may be driving the bulk of the star formation 
in these galaxies. This interpretation is corroborated by the large 
fraction of disturbed morphologies in the infrared luminous sam- 
ple, and also by our preliminary analysis of their local environ- 
ment: when compared to a matched sample of galaxies with no 
detectable infrared emission, infrared bright galaxies are more 
likely to have close neighbours, although the larger sample of 30 
LoCuSS clusters should be used to confirm this correlation. 
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